Enhanced angular overlap model for f-electron non-metallic systems 
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An efficient method of interpretation of the crystal field effect in non-metallic f-electron sys- 
tems, the enhanced angular overlap model (EAOM), is presented. The method is established on 
the ground of perturbation expansion of the effective Hamiltonian for localized electrons and first 
principles calculations related to available experimental data. Series of actinide compounds, An02, 
oxychalcogenides, AnOX and dichalcogenides UX2 where X = S, Se, Te and An — U, Np serve as 
a probe of effectiveness of the proposed method. An idea is to enhance the usual angular overlap 
model with ab initio calculations of those contributions to the crystal field potential, which cannot 
be represented by the usual angular overlap model (AOM) . The enhancement leads to an improved 
fitting and makes the approach intrinsically coherent. In addition, the ab initio calculations of 
the main, AOM-consistent part of the crystal field potential allows one to fix the material-specific 
relations for the EAOM parameters in the effective Hamiltonian. In consequence, the electronic 
structure interpretation based on EAOM can be extended to the systems of the lowest point sym- 
metries or/and deficient experimental data. Several examples illustrating the promising capabilities 
of EAOM are given. 

PACS numbers: 71.23.An 71.70.Ch 75.30.Gw 
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I. INTRODUCTION 



A list of phenomena traditionally discussed in the con- 
text of the crystal field (CF) like optical excitations, 
Schottky effect, Van Vleck susceptibility etc has been 
extended in last decades to include so intriguing and 
different manifestations of many-body effects as uncon- 
ventional superconductivity, Kondo-like behavior, mag- 
netic rearrangements or non-Fermi liquid. Complexity of 
these phenomena contrasts with apparent simplicity the 
one-electron, local CF potential in the effective Hamil- 
tonian. Nevertheless, its reliable determination meets 
serious difficulties even today, 75 years after Bcthc con- 
sidered the CF effect for the first timeii These concern 
not only calculations from first principles but also com- 
mon phcnomcnological schemes based exclusively on the 
symmetry arguments. In practice, only compounds of the 
highest symmetries can be handled satisfactorily. In most 
cases, the usual least squares fitting of all the parameters 
in the effective Hamiltonian requires additional, more or 
less heuristic reasoning referring to physical or chemical 
foundations. Alternatively, simplified phenomenological 
models are employed to circumvent the problem of over- 
numerous parameters. Accuracy of these models is one 
of the issues the present discussion addresses. We fo- 
cus on the angular overlap model (AOM) inspired by 
the molecular orbital theory^ and confirmed by fur- 
ther theoretical calculations!^ AOM assumes, similarly 
as the superposition model (SPM)£, the total CF poten- 
tial in the form of superposition of axial potentials due 
to ligands represented by a set of intrinsic parameters. 
The number of free parameters is remarkably reduced in 
comparison with the basal parametrization, especially for 
low-symmetry systems but the approximations are rather 
crude. A more refined analysis of the experimental data, 



having fundamental meaning for subtle magnetic proper- 
ties or cooperative phenomena at low temperatures, may 
require more rigorous methods. Can these simplified CF 
models be improved without increasing the number of 
free parameters? This question is discussed here on the 
grounds of ab initio calculations for a series of actinide 
compounds. 

There are two general ways to determine the CF poten- 
tial from first principles: one based on the perturbation 
theory and the idea of the effective Hamiltonian for lo- 
calized, open shell electrons^ and second developed on 
the grounds of the density functional theory (DFT)£& 
Taking into account proximity of the band states or/and 
widening of the localized nf states, DFT seems to be 
more suitable method for metallic systems. However, 
the early implementations of DFT based on the local 
spin density approximation (LSDA) had failed to pre- 
dict not only the subtle magnetic, low temperature ther- 
modynamic properties and the low energy spectra but 
even some of the main characteristics of the crystals. For 
instance, antiferromagnetic semiconductor UO2 becomes 
a ferromagnetic metal^ Over the years some inherent 
shortcomings of LSDA, like double counting of states or 
requirement for the electron density to be a slowly vary- 
ing function, were lessened or removed by various correc- 
tions: generalized gradient approximation, U approxima- 
tion, self-interaction correction or, more recently, the hy- 
brid density- functional theory^ Nevertheless, these im- 
provements, to our best knowledge, have not eliminated 
completely the problems related to the CF effect (see dis- 
cussion in Ref. 0, p. 206). 

An alternative approachjiSiii based on the perturba- 
tion expansion of the effective Hamiltonian seems to be 
more efficient and reliable in providing material-specific 
details of the electronic structure. The model was de- 
veloped for non-metallic systems successively by Sugano 
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and Schulmanji^ Newman and co worker s Faucher 
and Garcia?^ and others. (see Refs. Il0ll3l and references 
therein) The effective CF potential appears in this ap- 
proach as a sum of several contributions, most of which 
obey the assumptions of the mentioned simplified phe- 
nomenological models. 

In metallic systems, additional mechanisms, apart 
from those characteristic for insulating crystals, have to 
be considered^ mixing of the localized and band states 
(hybridization), static and dynamic screening of the con- 
duction electrons including an offset screening of the 
carriers occupying a virtual bond state. None of them 
can be regarded as a pure superposition of the nearest 
neighbor contributions. Moreover, the role of each mech- 
anism remains unclear; for instance, the hybridization 
term alone can represent almost the whole CF effect for 
some compounds^ but for others it is only one of the 
important contributions In addition, the dynamic cor- 
relations become increasingly important as the localized 
state nears the conduction band giving rise to the many- 
electron crystal field effects or even leading to a break- 
down of the effective Hamiltonian theory for localized 
electrons. Thus, an extension of the crystal field the- 
ory towards self-consistent models for mixed systems of 
localized and itinerant electrons attracts growing atten- 
tion recently^ However, material-specific results based 
on these efforts has not been achieved yet and we have to 
confine discussion to simpler systems in which an admix- 
ture of the band states to the localized ones or a mixing of 
the states localized on adjacent atoms can be treated as 
an additional perturbation term in our effective Hamil- 
tonian. We can only notice that there exists a subgroup 
of metallic systems for which the basic CF mechanisms 
remain essentially the same as those for ionic crystals^ 

Accuracy of the simplified models has been discussed 
on the grounds of numerical simulations for actinidc 
compounds^ It has turned out that the lattice contri- 
butions were not always negligible and the e$ contribu- 
tion (see section III Bp has not behaved as a characteris- 
tic, intrinsic parameter of the metal-ligand (ML) inter- 
action. Moreover, the present paper shows that the main 
AOM energies, e a and e^, may also fail a more rigorous 
test of their transferability. Thus we propose a quasi- 
phenomcnological approach that links the simplified phc- 
nomenological model with partial ab initio calculations: 
the enhanced angular overlap model (EAOM). 

The paper is organized as follows: the perturbation 
expansion leading to the effective Hamiltonian and the 
angular overlap model is outlined in the section [H] De- 
tails of the ab initio calculations, their results obtained 
for the first time for XJX 2 (X = S, Sc, Te)3£ and those 
reported previously for AnO^& and AnOX^L2L are dis- 
cussed in section. ITTT1 including their reliability, adequacy 
of the superposition approximation and variation of the 
intrinsic parameters with distance across the series. The 
enhanced angular overlap model is considered in sec- 
tion IIVI together with examples of its application. Con- 
cluding remarks are provided in section [3 To ensure 



self-consistency of the presentation some known formu- 
las used in the calculations are included in Appendices. 



II. FORMULATION 

A. Effective Hamiltonian 

Since the theoretical model applied here has been pre- 
sented elsewhere (see Ref. ^]and references therein), wc 
only recap its main points for clarity of further discus- 
sion. The basic assumption of the model is that all mag- 
netic electrons in a crystal occupy stationary orbitals, 5/ 
in the case of the discussed here actinide ions. The ini- 
tial infinite many-electron problem is reduced then to a 
single cluster consisting of a metal ion and nearest neigh- 
bors - ligands. The outside of the cluster is represented 
by the classic electrostatic potential. The cluster itself is 
treated as a system of weakly interacting groups of elec- 
trons localized on different ions. This allows one to ap- 
ply the group product function formalism and reduce the 
initial N-electron system to several subsystems of lower 
dimensions AifliSS The zero-order group product wave- 
functions are built up from the free-ion spin-orbitals ob- 
tained with the standard self-consistent Dirac-Slater pro- 
cedure and stabilization potential wall for anions, depth 
of which is determined by the Madelung energy. The 
function basis is restricted to the ground and the most 
important inter- and intra-ion excited electronic configu- 
rations. Projection of the initial function space onto the 
ground configuration subspace, contraction of the closed 
shell states, and renormalization due to the nonorthogo- 
nality of the wave-functions centered on different ions are 
the main steps in this approach. They lead to an effec- 
tive Hamiltonian defined in the restricted wave-function 
space spanned by the single 5/"-configuration states. 
The effective Hamiltonian contains several renormaliza- 
tion terms which can be regarded as a perturbation to 
the initial Hamiltonian. The theory may be easily veri- 
fied by the experimental data because the same function 
basis is employed in the conventional phenomenological 
description of the electronic structure (see Appendix A). 
A non-spherical part of the effective Hamiltonian defines 
the effective CF potential V(r) "seen" by a 5/-electron. 
Since all the wave-functions in our restricted function ba- 
sis have the same radial part -P(r), the Hamiltonian and 
all relevant operators can be contracted to the angular 
coordinates. The resulting operator of the CF potential 
V(r/r) is commonly written in a form of expansion in 

terms of the normalized spherical harmonics Cq ^ 

V{v/r) = Y,B kq Cf\v/r), (1) 

k.q 

where fc = 2, 4, 6 denotes rank of the spherical harmonic 
and q = — k, — k + 1, k runs over its components. B^ q 
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arc the integrals. 



B kg = / [-Pn/(r)] a V(r)C-W*( r /r)dr, (2) 



playing a role of adjustable CF parameters in the phe- 
nomenological theory. 



the right side of Eq. (J3J) with the matrix elements of the 
potential given by the expansion and using the prop- 
erties of tensor operators and the rotation matricesS^^ 
we obtain after some manipulation the following relation: 



B, 



kg 



WLe 



(6) 



B. Angular overlap model 

Angular overlap model (AOM) is a simplified phe- 
nomenological approach based on certain restrictive as- 
sumptions, inspired by the Hiickel molecular orbital 
model^^ From among various formulations appearing in 
the literature we chose one proposed by Schaffcr™ which 
does not refer to molecular orbital scheme. This for- 
mulation is consistent with the perturbation approach 
outlined above and the Newman superposition models 
discussed later. According to it. the CF potential V is a 
superposition of the independent contributions - poten- 
tials Vt generated by the nearest neighbors - ligands: 



7(r)~5^«t(r-R t ), 



(3) 



where Rt denotes the position of the t-th ligand. Addi- 
tionally, as an approximation, the local symmetry of the 
separated metal-ligand system is assumed to be axial. 

The AOM parameter ej, of the given ligand t is defined 
as a matrix element of the ligand potential vt in the co- 
ordinate system t with the z axis along the metal-ligand 
t linear ligator in which vt is diagonal: 



e* ft (iJ t ) = (± M K|±M), 



(4) 



where the index /i = 0(a), l(ir), 2(5), 3(^>) denotes the 
absolute value of the magnetic quantum number of the 
5/-electron. It is convenient to fix the energy scale by 



setting = 0. In practice we put = 



instead 



of 



Since only e M 's are used hereafter, we omit the 



tilde for convenience. Transformation properties under 
the rotation group -R3 of the I = 3 wave-functions allow 
one to express the matrix elements of V in the global 
coordinate system in terms of e^'s: 

(m |V| m>) = D !m (°. *t)DW„ (0, O f , $ t )4> 



(5) 

where D^m(0, 9t,$t) is the matrix element of the irre- 
ducible representation of the rotation group and 
Rt,@t, $i are the angular (global) coordinates of the lig- 
and t. Eq. JSJ) is the fundamental equation of AOM. It 
relates the matrix elements of the CF potential to the in- 
trinsic parameters describing the individual metal-ligand 
pairs through rotation matrices dependent on the geom- 
etry of the coordination polyhedron. 

For practical purposes it is advisable to relate the AOM 
parameters to the basic CF parameters B^q. Comparing 



where 



2k + 1 
7~ 



3 k 3 


3 /c 3 

—(JL /! 



(-l)^(2-^o) x 
£cf>*(O t ,$ t )4 (7) 



(8) 



and (:::) are the 3j symbols. The e M parameters are the 
mean values of the AOM parameters averaged over t. 
Their introduction into Eq. JHJ) is one of the possible so- 
lutions of the problem of several sets of intrinsic param- 
eters in the case of non-equivalent ligands. Note, that if 
W£ and given by Eqs. (|7I8|I are inserted into Eq. JfjJ), 
then the averaged e^'s cancel out. Moreover, the ratios 

can be treated as parameters of the model instead of 
e* 's. From the algebraic point of view, this is only a scal- 
ing of the parameters without loss of generality. In prac- 
tice, it is possible to estimate the ratios independently 
(this question is discussed later) and to consider the av- 
eraged quantities e M to be adjustable parameters. The 
W£ coefficients absorb all information about the geom- 
etry of the coordination polyhedron whereas the ratios 

encode differences in the AOM parameters due to the 
individual ML t distances. The distance dependence of 

has the exponential character—, yet, within a limited 
range of distances, it is fairly well approximated by the 
simple power function: 



(9) 



with the power exponents taking values in the range 
4.3 to 8.9 in the case of actinidc ions and simple 
Iigands£i2i2a 

An extension of the above equations to arbitrary num- 
ber of different anions is straightforward. Note that the 
partitioning © neglects the contribution from the out- 
side of the coordination cluster. This contribution and 
also other effects, which cannot be represented in the 
form of decomposition J3J), are included explicitly in the 
enhanced model presented below. 

AOM reduces the number of adjustable parameters de- 
scribing the CF effect to three. The strength of the model 
stems from the fact that the local inter-ion interaction 
parameters, e^, may be regarded as quantities charac- 
teristic for a given metal-ligand (ML t ) pair. They al- 
low one to verify unphysical solutions generated by false 
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minima of the fitting procedures in the standard para- 
metric analysis on the one hand and to indicate the en- 
suing approximations in the case of especially complex 
systems, on the other. Due to the dominant character 
of the renormalization terms (see section III CI and Eq. 
Ijl4|l ). the AOM parameters e a and e^ manifest several 
characteristic properties ii^ifliiiSliSi 

(i) their values reflect the spectrochemical ordering of 
the anions and 

(ii) decrease slightly with increasing atomic number 
along the lanthanide and actinide series and, indepen- 
dently, with decreasing oxidation number of the metal 
ion, 

(iii) e a > e„ > \eg\, 

(iv) e^/e* t const, for two different ML t and ML t ' 
systems (t ^ t'). 

(v) e^/e^. ~ const, for a given ML t pair, 

(vi) the e$ parameter, usually of minor importance, 
has been shown to be 'lattice sensitive— if obtained from 
the fitting of the experimental data. For this reason, 
transferability of es between various compounds seems 
to be questionable in the conventional AOM. 

Even though exceptions from the above rules happen 
(see the next section), they may serve as an instructive 
test of any set of CF parameters determined from the 
experimental data. This concerns not only the fitting 
results obtained with other approximate models but also 
the basic Bk q parameters, which can be "translated" to 
e n using Eq. 10 and a standard least squares procedure. 

As mentioned, one of the most widely applied approx- 
imate methods, the Newman superposition modclf£ii£ is 
based on the same assumptions as AOM. However, the 
role of intrinsic parameters is played by the Bk q param- 
eters for a separated linear ligator in a local coordinate 
system. Due to the assumed cylindrical symmetry of each 
ML subsystem, only the parameters with q = arc effec- 
tive. They are denoted hereafter as "bk" to distinguish 
them from the Bko parameters in the global coordinate 
system. The relation between the two sets of intrinsic 
parameters, AOM and SPM, can be easily obtained from 
Eqs. (|6I7I8[) by rewriting them for the specific case of the 
separated linear ligator: 



2k + 
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00 
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x E(-ir(2-^ )(4 o *) e ^ ( 10 ) 

The algebraic equivalence of the two sets of intrinsic pa- 
rameters evidenced by above relation allows one to ex- 
tend the literature data of their values, independently of 
the way in which they have been obtained. 

C. Contributions to the CF potential 

The procedure outlined in subsection lll Al lcads to sev- 
eral characteristic contributions to the CF potential V, 



corresponding to different mechanisms. Their discussion 
in next sections precedes a formulation of the enhanced 
angular overlap model. Thus, it is advisable to sepa- 
rate out these which obey the assumptions of the angu- 
lar overlap model (the AOM- consistent contributions), 
V AOM , and the residuum, V res : 

y _ yAOM + yres_ £qj 

In y AOM one can further distinguish the primary (V pr ) 
and renormalization (V ren ) components: 

yAOM _ ypr + yren_ ^ 

V pr represents the Coulomb interaction (direct and ex- 
change) of the ligand electrons, and the potential of the 
nuclei. It diverges from the point charge model (pcm)i 
due to the charge penetration (cp)2^ and inter-ionic ex- 
change (ex)^ effects included in V pr : 

ypr _ ypcm yep _|_ yex _|_ ypr.sh ^3) 

The intra-ionic excitations on the metal ion induced by 
the primary contribution lead to the AOM-consistent 
part of the shielding potential yv r - sh . 

yren com p r i ses t, ne m ain renormalization terms im- 
plied by the ligand-metal charge transfer excitations - 
the covalency contribution, V co ^ and non-orthogonality 
of the free-ion wave- functions localized on adjacent ions: 

the overlap contribution, V ov ^ and the contact shielding, 
ycsi±. 

yren _ yov + yes + yco ^ 

yres mc } uc j es a \\ ^ ne remaining terms. One can distin- 
guish y nn -P° l - the contribution of electric multipoles in- 
duced on ligands (polarization of nearest neighbors), V^ n 
- the electrostatic potential of the point charges and elec- 
tric multipoles induced on the all the ions outside the 
cluster (electrostatic potential of further neighbors) and 
yres.sh wn j cn symbolizes the shielding correction to these 
electrostatic potentials: 

yres _ ynn.pol _|_ yfn _|_ yres.sh (15) 

Due to the cancellation of various terms in the primary 
and residual contributions, V ren or V ov in essence, turns 
out to be the most important mechanism for the ionic 
compounds in favor of AOM^SiiSiiiiSI This corollary is 
supported also by the results presented in next section. 

Explicit formulas for all the above contributions are 
given in Appendix B. 

III. AB INITIO CALCULATIONS 

A. Details of the calculations 

The compounds under consideration represent a va- 
riety of the crystal structures. The point symmetry of 
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the metal ion varies from cubic one in UO2 (CaF2 struc- 
ture, space group Fm3m - 225), through D 2c i in VOX 
( PbFCl-type structure, PA/nmm - 129) C 2v in UTc 2 
(orthorhombic space group Immm - 71)|2i up to C s in 
/3-US 2 and /3-USe 2 (PbCl 2 -type structure, Pnma - 62) M 

Zero-order free-ion wave-functions have been generated 
using self-consistent Dirac-Slater procedure ATOMi2i2i 
with the stabilizing potential well determined by the 
Madclung energy^ The calculations have been per- 
formed in the crystallographic coordinate system. The 
lattice sums of static and induced multipoles generated 
by the set of the crystal electrostatic equilibrium equa- 
tions l|B13jl in Appendix 151 have been calculated accord- 
ing to Eq. (|B12|) using a modified version of the pro- 
gram CHLOEii^iSi and multipole polarizabilities from 
Ref. The summation in the multipole expansion 

(|B12f) was limited to monopoles, dipoles and quadrupoles 
(n = 0,1,2 ) in the present calculations. The effect 
of the outer electrons occupying the 6s and 6p closed 
shells of the metal ion have been estimated via Stern- 
heimer's shielding factors, tjfc, scaling the corresponding 
Bkq , a^^. The charge penetration, exchange and renor- 
malization terms have been calculated using the program 
LF developed for the /-electron compounds^ 

The free-ion part of the effective Hamiltonian contains 
the intraionic Coulomb repulsion of the 5/ electrons of 
the U(4+) ion, controlled by the Slater integrals F k , 
k = 2,4,6, spin-orbit coupling with the Cs/ parameter 
and further corrections of higher order (see Appendix 
lA"|) . The "free-ion" parameters depend on the crystal, in 
which the metal ion is embedded, and may vary in certain 
range, modifying the inter-term spacing^! Fortunately, 
these differences are limited and they have rather neg- 
ligible impact on the splitting of the ground multiplct. 
Thus, we adopt the values obtained earlier from the fit- 
ting of the optical spectra for U 4+ in ThGe04 (in meV)4£ 
F 2 = 5339, F 4 = 4833, F 6 = 3018, a = 3.72, f3 = -81.8, 
7 = 148.8, M° = 0.124, M 2 = 0.069, M 4 = 0.047, 
P 2 = 62.0, P 4 = 31.0, P 6 = 6.2, Cs/ = 224.2. The calcu- 
lations of the eigenvalues and eigenvectors and/or fitting 
of the measured energies of electronic transitions have 
been performed with a set of f-shell empirical computer 
programs developed by M. F. Reid^ supplemented with 
a subprogram POD 4 ^ for calculation of some thermody- 
namic quantities. 



B. Reliability of the results 

This type of ab initio calculations was performed pre- 
viously for, among the others, UO2 and UOS ^L2L^ the 
compounds for which comprehensive experimental data 
are available, including optical and/or inelastic neutron 
scattering (INS) spectrai 44 i 45 i 46 Discrepancy of the ob- 
tained theoretical Bk q parameters and the phenomeno- 
logical values determined by a fitting to the experimental 
data did not exceed 20-30%. Now the calculations are 
extended to the UA2 sub-scries. For these compounds 



the published experimental data arc not so richi^aLaSi 
We can reproduce certain thermodynamic averages, e.g. 
the paramagnetic Van Vleck susceptibility (see Appendix 
EJ) depicted in Fig. ^ The corresponding energy levels 
are shown in Fig. |2 and they are discussed later. The 
experimental slope, shape and the variation of the tem- 
perature independent susceptibility at the lowest temper- 
atures along the series agree quite well with the results 
of the ab initio calculations. A discrepancy between the 
theoretical and experimental lines visible for USe2 and 
UTe2 can be attributed to the intcrionic exchange not 
included in the model calculations. Judging from the 
mutual shifts of the lines, the interaction has to be of 
ferromagnetic type for USe2 and greater in the absolute 
value than the antifcrromagnctic type observed for UTe2. 
The first excited level lying at 8.0, 1.6 and 2.5 meV for 
U 4+ ion in US2 , USe2 and UTe2 , respectively (see Fig. |5J) 
influences the magnetic properties of these systems with 
a singlet ground state. Above 20 K, UTe2 and especially 
USe2 behave as if they had a degenerate ground state. 
Relatively weak U-U exchange interaction in USC2 (U-U 
distance of 0.423 nmS) turns out to be sufficient to in- 
duce the long-range ferromagnetic order below 14 K4£ 
Note, that the first excited state has the lowest energy 
just for this compound. USe2 is the only dichalcogenide 
that orders magnetically^ No ordering in UTe2 with the 
shortest U-U distance 0.378 nm may be due to the fact 
that only one of seven neighboring uranium ions is placed 
at that distance, whereas the remaining four are placed 
at 0.490 nm and two at 0.416 nm. A different curvature 
of the theoretical and experimental reciprocal suscepti- 
bilities observed for USe2 below 60 K may be ascribed 
to the magnetic fluctuations increasing with temperature 
approaching the critical point, which were not taken into 
account in the model calculations. The effective magnetic 
moments, of 2.94 BM, 3.01 BM and 3.09 BM, are lower 
than those derived from the experimental curves, of 3.25 
BM 3.20 BM and 3.21 BM, reported for US 2 , USe^ and 
UTc2r^ respectively. 

Figure [5] shows the splitting of the ground term 3 H^ of 
the t/ 4+ (5/ 2 ) ion in the UX2 series, obtained by simulta- 
neous diagonalization of Hamiltonian (|A1IA2|) , with the 
CF parameters determined from first principles shown in 
Table [I] and the free-ion parameters listed in Appendix 
lAl These results have been employed in calculations of 
the temperature dependencies of the magnetic suscepti- 
bilities discussed above. Bk q parameters in Table U reveal 
no regularity along the series that might be expected, for 
instance, from the spectrochemical ordering of ligands. 
This observation is not only a consequence of the inher- 
ent ambiguity of the CF paramctrization in low symme- 
try systems^ It would be difficult to find any trend also 
for the energy levels shown in FigEI Analogous behav- 
ior of the AnOX series has been ascribed to an influ- 
ence of the competing oxygen and chalcogenide groups 
in the W£ coefficients Q.-SLSI Now it becomes evident 
that also in the case of one type of anions the coordina- 
tion geometry may obscure the expected regularities so 
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FIG. 2: Splitting of the uranium (4+) ground term 3 H4 in 
the XJX2 crystals obtained by simultaneous diagonalization 
of Hamiltonian l|AllA2fl with the CF parameters determined 
from first principles. The lower part of the diagram uses the 
scale enlarged several times as compared with that for the 
upper part. 



clearly manifested, on the other hand, by the intrinsic 
AOM parameters (cf. Table HT|) . As we see in next sec- 
tion, possible trends in CFP sets revealed elsewhere using 
an independent method based on certain conserved quan- 
tities associated with CF parameters^ may be governed 
by the AOM-consistent part of crystal field. 



C. AOM-consistent part of crystal field 

The results displayed in Table [i] are indicative of an 
essential meaning of the AOM-consistent contributions. 
Nevertheless they also show that the residual contribu- 
tion may become crucial for some parameters as in the 
case of B44 for US2 or B22 for UTC2 . More detailed data, 
which are presented in Table [I] for USe2, give some idea 
about the role of the particular CF mechanisms. Mutual 
compensation of the various primary components and im- 
portance of the renormalization terms is clearly mani- 
fested there. Handling each non-equivalent ML pair 
independently in compounds like US2 (with six different 
ML t distances and two non-equivalent crystallographic 
positions of the sulphur ion) would multiply the number 
of the intrinsic parameters, making the model practically 
intractable. Usually, the dependence of the intrinsic pa- 
rameters (or their renormalized counterparts s M ©) 
on the ML distance, R, is assumed to have a definite 
characterful. In a limited range of distances around an 
average value for a given ML bond, it is approximated 



by the power function (j5J). The power exponents, a^, are 
treated then as an additional characteristic of the ML 
bond that allows one to maintain the minimal number of 
independent parameters. 

The simulations of the (R) functions have been per- 
formed for all the uranium linear ligators occurring in 
the series U 4+ - A 2 ". As Eqs. (|B4IB10j) in Appendix El 
show, the functions are determined by squares of metal- 
ligand overlap integrals and Madclung energies of the ions 
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TABLE I: AOM-consistent \V2l and residual 115t contributions to the CF parameters calculated from first principles for the 
UX2 series in the crystallographic coordinate systems^- The primary contribution V pr l|13[l . renormalization term V ren (SJ, 
ligand polarization y nn -P° l 1B12I I and potential of the further neighbors V^ n are specified for USe2 in parentheses. All values 
in meV. 
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-3 


-206 


-200 


(-194 - 5) 


-6 


(-1 - 5) 


-206 








B22 


-79 


-6 


-85 


-24 


(-26 + 3) 


-3 


(-9 + 6) 


-26 


28 


42 


70 


B40 


-88 


2 


-86 


-71 


(-52 - 19) 


4 


(4-0) 


-67 


-436 


-184 


-619 


B41 


82 


6 


88 


39 


(32 + 7) 


-5 


(0-5) 


34 








B42 


-274 


-8 


-282 


-264 


(-297 + 33) 


-16 


(6 - 22) 


-280 


364 


79 


442 


B43 


257 


45 


301 


186 


(236 - 49) 


26 


(25 + 1) 


212 








B44 


-60 


-51 


-111 


-64 


(-62 - 2) 


-48 


(-21 - 27) 


-112 


74 


-26 


48 


Bm> 


174 


-2 


172 


78 


(80 - 3) 


1 


(2-1) 


78 


87 


4 


91 


B 6 i 


-233 


-1 


-233 


-230 


(-311 + 82) 


-1 


(-1 + 0) 


-231 








B(,2 


394 


5 


398 


340 


(441 - 101) 


3 


(-5 + 9) 


344 


11 


-6 


6 


B&3 


171 


-1 


170 


203 


(248 - 45) 


-1 


(-1 - 0) 


202 








Bq4 


37 


18 


55 


34 


(59 - 26) 


9 


(8 + 1) 


43 


-122 


26 


-96 


Bqt, 


-369 


-10 


-379 


-320 


(-383 + 63) 


-4 


(-4 - 0) 


-324 








Bee 


243 


9 


252 


212 


(262 - 50) 


6 


(2 + 4) 


219 


-151 





-151 



TABLE II: Results of ab initio calculations of the AOM pa- 
rameters (in meV) corresponding to the AOM-consistent part 
of the CF potential for various An 4+ — X 2 ~ systems and the 
average inter-ion distances R av . 





Rav [nm] 


e CT 




es 


xj i+ -6 2 - 


in: 








uo 2 


0.237 


350 


209 


64 


uos 6 


0.234 


340 


211 


63 


UOSe 6 


0.236 


324 


196 


55 


UOTe 6 


0.237 


316 


186 


46 


u 4 +-s 2 - 


in: 








uos 


0.293 


211 


96 


32 


us 2 


0.289 


265 


102 


36 


U 4+ -Se 2 ~ 


in: 








UOSe 


0.304 


204 


91 


31 


USe 2 


0.301 


240 


88 


30 


Tj4+_ Te 2- 


in: 








UOTe 


0.325 


185 


90 


31 


UTe 2 


0.317 


190 


88 


29 


Np 4+ -0 2 " 


in: 








NpOa 


0.235 


317 


186 


56 


NpOS 


0.233 


304 


181 


53 


NpOSe 


0.235 


291 


169 


47 


Np 4 +-S 2 " 


in: 








NpOS 


0.291 


191 


87 


30 


Np 4 +-Se 2 " 


in: 








NpOSe 


0.302 


183 


82 


29 



b From Ref. El 



in crystals. To ensure comparability of the results ob- 
tained for different ML pairs, a virtual ML2 crystal of 
the CaF 2 structure has been employed, where ligands 
form a cube with its center occupied by the metal ion. 
The Madclung potential at the L and M sites in this 
structure is given by the formulae: U L = — 8.14e/a and 



U = 15.13e/a, respectively, whereas the lattice con- 
stant a is related to the ML distance R by the expres- 
sion R = \/3a/4:. The intrinsic parameters depend on 
U L and U M in a non-trivial way through the zero-order 
wave-functions and explicitly due to the renormalization 
terms <|Ij4l IB6II . This implies the use of different free- ion 
wave-functions for each R. The calculations have been 
performed for seven values of R distributed uniformly 
between 0.18 nm and 0.40 nm. The radius D of the sta- 
bilizing potential well for the negative ions of the U L 
depth has been assumed to be equal to —2e/U L . 

The results are shown in Fig. [3] Due to the predom- 
inating renormalization terms (cf. Fig. the e^R) 
functions have an exponential character for the ML dis- 
tances around the average values. For larger distances, 
where the electrostatic contributions dominate, they take 
a form of a polynomial. The slope of e CT (R) decreases in 
the limit of short distances, more visibly on going from 
the oxide to tclluridc. The function e a (R) reaches even 
a maximum for the U 4+ -Tc 2_ bond. Fig. 0] indicates 
that the importance of the charge penetration rapidly 
increases as compared with that for other mechanisms. 
Since this contribution has opposite sing for e a and e T , 
the ratio e^l 'e a increases in the limit of small distances. 

A similar behavior, although not so evident, can be 
deduced from the data presented in Table II of Ref. 
for the Pr 3+ -Cl _ , if the Coulomb contributions arc scaled 
using the Stcrnheimcr shielding factors and the SPM pa- 
rameters are converted into the AOM ones by Eq. H10[) . 
An open question is whether the observed increase of 
the e n /e a ratio at the lowest distances characteristic for 
dense systems is a true property of metals. We may only 
note that analogous distance dependencies of the AOM 
parameters for the Sm 2+ -Cl~ and Sm 2+ -F~ systems de- 
rived from the data reported in Ref. did not confirm 
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FIG. 3: Dependence of the inter-ion effective interaction pa- 
rameters of the angular overlap model on metal-ligand dis- 
tance for the U 4+ — X 2 ~ pairs in a hypothetical crystal of 
CaF2 structure. The circles, squares and triangles represent 
actual intrinsic parameters e a , and e$ obtained for UX2. 

this conjecture. 

The model calculations taking into account the varia- 
tion of Madclung energy with the ML distance lead to 
lower values of the power exponents as compared with 
earlier estimations^ The AOM parameters determined 
for the actual crystals AnOX and AnX2 are listed in 
Table ITT1 and shown in Fig. [3J Each ML distance occur- 
ring in these compounds is represented in Fig. |3|by three 
points corresponding to the e a , e^, and e$ parameters. 
We see that they do not follow the smooth lines discussed 
above. The crucial e a parameter turns out to be simul- 
taneously the most irregular one. This is because it re- 
flects the discontinuity of the Madclung potential most. 
In light of the present calculations, the properties (i-v) of 
the AOM parameters specified in section III Bl should be 
treated with caution not only because of very existence 
of the residual contributions but also because of the re- 
vealed here inherent irregularity of the AOM-consistent 
contributions. 



IV. ENHANCED ANGULAR OVERLAP MODEL 

A. Formulation 

In view of apparent irregularity of the AOM parame- 
ters and the erratic behavior of the residual contribution 



(|15fl discussed in the preceding section and manifested in 
Fig. |2]and Tablc[I]we propose to link the phenomenolog- 
ical AOM approach with the ab initio calculations in the 
enhanced angular overlap model (EAOM). Our model 
assumes each ordinary CF parameter to be composed of 
two components: the main, adjustable one, parametrized 
according to Eqs. (Itil.Sti and the fixed residuum, cor- 
responding to the potential (|15fl . determined from Eq. 
l|Bl2)D : 

B kq = Y. W k^ + B l7 ( 16 ) 

with playing the role of the phenomenological EAOM 
parameters. The separation of the residual, off-AOM 
contribution from the global CF parameters supports 
a more precise and consistent AOM paramctrization of 
the remainder. In consequence, the intrinsic charac- 
ter of the corresponding EAOM parameters is genuinely 
maintained. EAOM is not much more complicated than 
the parental AOM but the main inherent shortcoming 
of the latter is removed by excluding explicitly the off- 
AOM contributions from the simplified phenomenolog- 
ical treatment. These off-AOM contributions are esti- 
mated in our model from the first principles. Note, that 
very construction makes the model exact (at least in the 
frames of the one-electron approximation) provided the 
off-AOM contributions are determined precisely, what, of 
course, is hardly possible. 

The AOM-consistent contributions listed in Table 
ITT1 can be regarded as a crude theoretical estima- 
tion of the EAOM parameters. Their more accurate 
determination requires an involved self-consistent ap- 
proach to account for the energy-dependent parameter- 
renormalization terms^ Thus, treating this part of the 
CF potential as a phenomenological quantity allows one 
to handle these self-energy effects in the simplest possible 
way. 

Naturally, the EAOM parameters obtained from the 
fitting of the observed electronic energy levels cannot 
have the same values as the AOM parameters derrived 
from the analogous fitting within the conventional AOM 
approach. Moreover, the rules (i-v) observed for the 
usual AOM parameters (see section Hi Bll seem more jus- 
tifiable in the case of the theoretical, AOM-consistent 
contributions and thus, also for the EAOM parameters. 
Nevertheless, in light of the calculations presented in sec- 
tion llll CI still they remain acceptable in a rather limited 
range of the ML distances, in vicinity of their averaged 
values. 

Compounds containing several groups of symmetry- 
equivalent ligands require the ratios in the geomet- 
rical coefficients W£ q J7J) to be determined. The ob- 
served above (see section IIII CI and Fig. irregularity 
of the EAOM parameters due to, among the others, the 
Madelung energy, can be reproduced in the model by 
employing the ab initio values of these ratios. Thus, it is 
advisable to calculate s^'s for each individual ML t pair 
rather than applying the approximation 
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FIG. 4: Distance dependence of the main contributions to the e M parameters: point charges (pcm), overlap (ov), covalency 
(co), charge penetration (cp), exchange (ex), contact shielding (cs). 



The same theoretical calculations may also be 
employed to formulate further. more restrictive 
paramctrizations, which may be helpful in the case 
of the most complicated experimental data. For 
instance, fixing both the ratios 



(17) 



equal to their theoretical values leads to a single- 
parameter version of the model: 



B kq = BH S + W kq e 



(18) 



with 



Wi 



kq 



2k + 1 



3 k 3 


3 k 3 
-/i /i 



£(-l)"(2-^ ) x 

X;^ fe) *(e t ,$ t ) s >^(i9) 



The e parameter in Eq. (|18f) corresponds directly to 
e CT but also the and contributions enter into the 
model through the s^ ia ratios in the above Wk q coeffi- 
cients. Hence, a variation of e modifies also the e n and e$ 
contributions to the crystal field in proportions given by 
s^o-. Note, that Eq. (|18f) without the residual part would 
be a simple scaling of the crystal field effect predicted by 
the AOM-consistent part of the ab initio calculations. 

An intermediate, two-parameter version of the model 
can be defined in several ways by fixing any of the s M<T 
ratios or their combinations. These single- and two- 
parameter versions, should not be confused with mod- 
els commonly employed for the transition elements^^ 
which omit simply the it and S contributions or only the 
d contribution. 

Generalization of the model and the equations (|16I19|) 
to systems containing different ligands is straightforward 
(see the example given below) . 



B. Applications 

Due to the highly reduced number of the free param- 
eters (up to the single-parameter version, Eq. H18fl ). the 
simplified EAOM parametrizations of the CF effect may 
be especially helpful in all these cases of incomplete or 
more complicated experimental data so frequently met in 
the most interesting /-electron systems. Three possible 
types of applications of EAOM are exemplified in what 
follows. 

The first one, already discussed in Ref. [^J concerns 
interpretation of certain magnetic properties of NpOA 
where X = O, S and Se. The procedure might be seen 
as a general method of interpretation of the electronic 
properties based on transferability of the intrinsic pa- 
rameters between different compounds, using the phe- 
nomenological CF parameters for the specimens assumed 
to be known, i.e. the compounds, the electronic struc- 
ture of which is believed to be reliably determined. In 
the example under consideration, the parameters (Bk q ) u 
reported for UO2, VOX, the ones of the most widely in- 
vestigated uranium compounds, were employed^2ii4 The 
corresponding AOM parameters e^f and e^ x were de- 
termined from system of equations (|16f) and its version 
adapted to two different ligands, oxygen and chalco- 
genide, 



Bkq — B 



kq 



J V kq e M 



kq e M 



(20) 



exploiting the ab initio values of the ratios. In the 
second phase, the procedure was reversed to estimate un- 
known CFP's for the less explored experimentally nep- 
tunium oxychalcogenides. First, the e^ p and e^ pX pa- 
rameters were derived by scaling their uranium counter- 
parts in terms of the ratios of the corresponding the- 



oretical values: (er[ p 



The B^ p parameters 
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obtained, again, from Eqs. I|16l I20f) were, in the final 
step, adjusted to match the experimental values of the 
ordered magnetic moments. In this phase, a simultane- 
ous diagonalization included, apart from the CF Hamil- 
tonian and the "free-ion" interactions, the intermetallic 
exchange interaction in the zero-temperature mean-field 
approximation. The magnetic properties of the neptu- 
nium oxychalcogcnidcs, including the anisotropic ground 
state magnetic moment, the temperature dependencies of 
the paramagnetic susceptibility and the magnetization at 
K, were described satisfactorily in this approach (see 
Ref. |^for further details). The example illustrates effi- 
ciency of the method based on EAOM in interpretation 
of the electronic structure in the case of inconclusive ex- 
perimental data. The intrinsic parameters collected in 
Table ITT1 may serve as a tentative source of the relations 
between them in further applications of EAOM. 

The second example concerns interpretation of the 
thermodynamic and magnetic properties' 4 ^ and inelas- 
tic neutron scattering (INS) spectra of UOSi 46,56 In the 
analysis of the experimental data^Sl the initial energy 
levels assignment was deduced on the grounds of first- 
principles calculations. The fitting of the INS data^^i 
was performed in two steps. First, the AOM was applied 
to obtain the identified INS transition energies. Sub- 
sequent refinement of the corresponding usual CF pa- 
rameters yielded agreement with the remaining experi- 
mental datai 4 ^*' 4 *" A special attention was paid to the 
ordered magnetic moment fi or d and relative intensities 
of the INS transitions. However, the problem with the 
infinite number of acceptable solutions could not be re- 
solved. Namely, the experimental data for this tetrag- 
onal Ci v system could be satisfactory reproduced with 
any value of the B20 in the range of —50 to —223 meV, 
provided the remaining parameters, -B40, B44, Bqq, Bq4 
were simultaneously adjusted to the given i?2o- Some 
features of the electronic structure were varying with the 
CF parameters in this ambiguity range but they were not 
detectable in fi or d and in the position and shape of the 
most intense INS lines below 100 meV. In addition, it 
was demonstrated 2 ^ that the INS excitations to two of 
the highest levels of the 3 i?4 term, T 2 and r 1 ^ 2 ), could 
be invisible if B20 had decreased below -146 meV under 
actual incident neutron energy and angle, just as it was 
observed. Eventually, the set of B kq parameters corre- 
sponding to the threshold value of B20 = — 146 meV was 
accepted as the closest one to the initial estimation. 

The above problem disappears if EAOM is applied 
with its inherent constrains. To start with, the single- 
parameter version of EAOM, Eq. (16), is adapted to 
account for the two different anions occurring in the co- 
ordination sphere - the oxygen and sulphur: 

B kq = B^ + W°e° + W k s q e s (21) 

with coefficients W given by Eq l(T5)l. B T k & q and the re- 
maining data taken from Ref. l20t It turns out that this 
simplest version of EAOM, with only two effective pa- 
rameters e and e , describes fairly well the three elec- 



tronic INS transitions^^ 74 meV (r 5 (l) -> T 3 ), 83 
meV (r 5 (l) -► T 4 ), and 87 meV (r 5 (l) -> T 5 (2)). These 
intervals can be easily matched by a minor refining of 
the ab inito ratios s WCT . Figure |3] (c) shows the result- 
ing energies and the transition probabilities estimated in 
terms of the squares of the matrix elements of the Zee- 
man operator between the initial and final states^ The 
Debye- Waller and 5/-electron form factors are not taken 
into account since we are interested only in comparison of 
various numerical simulations without a direct reference 
to the experimental recordings. The Gaussian shape of 
the figured transitions have only an illustrative charac- 
ter. The present simulation of the INS transitions does 
not differ much from the previous one shown in Fig. [S] 
(b)i2£ Taking into account a finite width of the experi- 
mental lines and limited range of the measurable ener- 
gies both the simulations seem acceptable. In particu- 
lar, the main lines around 80 meV are quite similar to 
the measured onesi 4 ^^ The line at about 40 meV was 
neither observed nor excluded by the experimental data 
since it is relatively weak and lies in the region of the 
phonon sidebands. The ordered magnetic moment of 2.11 
MB as determined from the obtained wave- functions, is 
higher than the observed one of 2.00 MBiS Neverthe- 
less it can be reduced to about 2.04 BM due to the over- 
lap and covalency effects^ Our phenomenological results 
can be compared with ab inito predictions. The param- 
eter e° = 370.2 meV is about 10% higher than the theo- 
retical value of e CT from TablclTTl whereas e s = 157.1 meV 
is lower than its counterpart. The differences between the 
corresponding phenomenological and theoretical parame- 
ters become even more pronounced in the case of the B kq 
parameters listed in Table lTTTl Taking into account inher- 
ent limitations of the ab inito calculations due to numer- 
ous approximations and uncertainty of the Sternheimer 
shielding factors, dipole and quadrupole polarizabilities, 
such a divergence seems inevitable. We assume the ratios 
of the theoretical values of the EAOM parameters, J§J 
II 711 , used in formulation of the approximate phenomeno- 
logical models, Eqs. (|16l I18f) . to be more reliable than 
the theoretical parameters themselves because of an ex- 
pected partial cancellation of the calculation errors. The 
present EAOM values of the B kq parameters lie between 
those obtained using the conventional AOM and the re- 
fined model from Ref. |2(1 The divergence between the 
sets of the phenomenological parameters, especially be- 
tween these in the two last lines of Table ITTT1 is moderate. 
Further INS investigations with higher incident neutron 
energies could decisively verify the predicted positions of 
the T 2 and T 1 ^ levels. 

The low symmetry crystals, UA2, considered in sec- 
tion IIII Bl may serve as a next example of application of 
EAOM. The model allows one to reduce 15 (US 2 , USe 2 ) 
or 9 (UTe2) B kq parameters specific for the point group 
symmetries in these compounds to only 1-3 EAOM pa- 
rameters. Note, that very construction of any version of 
EAOM, including the single-parameter one, Eq. i|18|) . en- 
sures accuracy of the electronic structure simulation not 
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FIG. 5: Simulation of the INS spectra of UOS in various 
crystal field models: (a) - ab initio calculations, (b) - refined 
model from Ref. H3, (c) - EAOM. 



TABLE III: Comparison of various Bu q sets obtained for UOS 
(in meV). 





B20 


B40 


B44 


Beo 


£>64 


Ab initio calculations 


-180 


-545 


-119 


424 


104 


AOM 


-117 


-567 


-11 


632 


250 


Model from Ref.20 


-147 


-579 


-39 


595 


447 


EAOM 


-159 


-624 


-21 


627 


380 



worse than that obtained from the ab initio calculations. 
This is because the ab initio calculations determine not 
only the starting values of the parameters but also the co- 
efficients W and residual off- AOM terms in Eqs. (|16I18|I . 
Thus, with EAOM one can try to reproduce the elec- 
tronic structure of any compound even if the available 
experimental data arc limited merely to a single ther- 
modynamic characteristic like the magnetic susceptibility 
discussed in section llll Bl In the specific case of VX2 , the 
theoretical curves displayed in Fig. ^ represent simulta- 
neously the initial phcnomcnological dependencies in our 
model. The corresponding EAOM parameters listed in 
Table [H] are thus the natural and right starting values 
in subsequent steps of the further interpretation of the 
electronic structure. 



V. CONCLUDING REMARKS 

The crystal field potential in the quasi- 
phenomenological enhanced angular overlap model 
(EAOM) proposed here is divided into two parts: the 
main part, adjustable, which comprises the contributions 
amenable to the conventional AOM paramctrization and 
the residual one, fixed, representing all the remaining 
terms. The fixed, off-AOM residual part contains 
the electrostatic contributions of further neighbors 
and the ligand polarization contribution, accuracy of 
determination of which hinge on the available ionic 
polarizabilities and the shielding factors. Much more 
complicated renormalization terms in the CF potential 
are handled in the phenomenological way together with 
the core, charge penetration and interionic exchange 
contributions. The paramctrization applied is known 
from the usual AOM approach. Its simplicity is a 
consequence of the axial local symmetry of the encoded 
contributions and their additivity. 

Aptness of such a discrimination between the contribu- 
tions is illustrated by an increasing amount of examples 
of first principles calculations in the literature and also 
by the results presented in this paper. Additionally, the 
calculations provide the ratios of the intrinsic parameters 
which may be employed in further applications of EAOM. 
The simulations performed for the virtual XJX2 crystal of 
CaF2 symmetry give an idea about dependence of the 
EAOM parameters on metal-ligand distance. They show 
the evolution of the mutual relations between the AOM- 
consistcnt contributions. The observed behavior of the 
AOM-consistent contributions to the intrinsic parameters 
in the limit of the shortest metal-ligand distances touches 
a more complex problem of the CF effect in metallic sys- 
tems. A comparison of these simulation with the ab ini- 
tio calculations for the actual VX2 compounds evidences 
certain volatility of the intrinsic parameters. It turns 
out that the assumption, commonly accepted for models 
based on the superposition anzatz (conventional AOM or 
SPM and their various modifications), saying that the in- 
trinsic parameters are smooth functions of metal-ligand 
distance is not always valid due to various Madelung en- 
ergies of the ions in the real crystals and a high sensitivity 
of the intrinsic parameters to these energies. This points 
to necessity of enforcing each particular implementation 
of the phenomenological model with the ab initio calcu- 
lations for an actual specimen just as it is postulated in 
EAOM. 

EAOM ensures a highly compact description of the 
electronic structure of even so complex ionic systems as 
the actinide crystals. The examples discussed in the pa- 
per represent various kinds of conceivable applications. 
The case of experimental data limited to only certain 
magnetic characteristics of systems of low crystal sym- 
metries is illustrated for the UX2 series. In other exam- 
ple, the electronic structure for the NpOX is predicted by 
transforming the intrinsic parameters for the correspond- 
ing VOX compounds and exploiting relations between 
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the parameters calculated ab initio. Various experimen- 
tal data are available in the case of UOS. They include 
the magnetic and thermodynamic properties as well as 
the inelastic neutron scattering spectra. Nevertheless, 
the electronic structure cannot be resolved within the 
conventional CF parametrization scheme without addi- 
tional assumptions of rather heuristic nature. EAOM is 
shown to be capable of describing of all of the observed 
properties unequivocally. 

EAOM is open for further improvements. We mention 
here about a more accurate potential that could be ob- 
tained from the lattice self-consistent calculation; 



to 

correct the zero-order wave- functions used in the calcula- 
tions. Direct evaluation of multipole polarizabilities and 
shielding factors for a given specimen may be important 
for calculations of the residual part of the crystal field, 
since these quantities depend on the Madelung energies of 
the ionsi2£ Perhaps, an analogous idea of partitioning of 
the CF potential into phenomenological and fixed parts 
may be also applied to much more complicated metallic 
systems. 
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APPENDIX A: 



(3 and 7 parameters are associated with the two-body cor- 
rection terms. G(G2) and G(Rr) are Casimir operators 
for the G2 and R7 groups and L is the total orbital an- 
gular momentum. The electrostatically correlated spin- 
orbit perturbation is represented by the Pk parameters 
and those of the spin-spin and spin-othcr-orbit rclativistic 
corrections by the Mj parameters. The operators associ- 
ated with these parameters are designated by rrij and pk 
respectively. 



APPENDIX B: 

Main contributions to the one-electron CF potential 
in non-metallic crystals are listed here. For derivations 
and a detailed discussion see Ref. [ToL Some numerical 
questions related to multicenter integrals and summation 
of weakly convergent infinite series arc dealt with in Rcfs. 
ITllll2a 

From definition, the AOM-consistent part of the CF 
potential and thus the AOM-consistent contributions in 
Eqs. 113I14|I are partitioned due to ligands according to 
Eq. (J3J. Therefore, without lose of generality we can 
consider only single ligand potential vt in the local coor- 
dinates centered at the metal site, the z axis of which is 
directed towards the ligand. Due to the axial symmetry 
of the ligand potential the orientation of the x and y axes 
is immaterial. Atomic units arc applied throughout this 
section. 



Rt 



v?{v) = E 



J(Xtr(r),Xtr(r)) 



Ir - Rt 



(Bl) 



(B2) 



The parametric Hamiltoniar*i2i2&3i contains the free- 
ion, spherically symmetric part, H$, and crystal-field po- 
tential, V, 



(Al) 



where the summation index i runs over all /-electrons. 
The free-ion part can be written as follows: 



H 



E ave + E Fk {nf,nf)f k + ( 5f A 30 

k=2,4,6 

aL(L + 1) + (3G(G 2 ) + jG(R 7 ) + 
E M 3 m 3 + E /"'• 



Pk 



(A2) 



k=0,2,4 



k=2,4,6 



where E ave is the spherically symmetric one-electron part 
of the Hamiltonian, F k (nf,nf) and £5/ represent the ra- 
dial integrals due to the electrostatic and spin-orbit in- 
teractions, while fk and A$o are the angular operators 
corresponding to these interactions, respectively. The a, 



v? 



(r) = -^AT( XtT (r), XtT (r)) 



(B3) 



C = Y\^r){Xtr\[^\ho\ V )-2h + 

T 

+ £(^--%>V))IXfT'><Xfr'fl (B4) 

t>,T> 

E Eti&ixtr) i 2 [2 6) - m v , &,)] - 

V T 



t',r' 



x [2J(xtT,Xt>T>) - K(xtr,Xt'r')]} 



E 



htr\Xtr){Xtr\hu 



(B5) 



(B6) 



13 



ip m and £ v in above equations denote the 5/ and 6s, 6p 
orbitals of the metal ion obtained from the Dirac-Slater 
calculations for the free ion, xtr stand for the analogous 
ns and np orbitals of the anion and m, is, t are the 
corresponding sets of quantum numbers. J and K are 
the usual Coulomb direct and exchange operators (see 
for instance Ref. ITflh . 



vf cm (r) represents only ligands as point charges 
(monopoles) . 

The remaining contributions are convenient to be pre- 
sented as a part of the electrostatic model, which includes 
all ions in the crystal represented by a sequence of point 
monopole, dipole, quadrupole etc.: 



ho 



yfr. 



(B7) 



where v^f and v^ t are the Dirac-Slater free-ion metal and 
ligand potentials. 

h tT = h tT -J2 Yl \Xt'T>)(Xt'T>\ht T \Xt»T»)(Xt"\ (B8) 

t'r' V't" 



h tT = ho- J(xtT,Xtr) (B9) 



A tT = (<p\htr\<p) - (XtMXtr) (BIO) 



V 



ynn.pol , yfn , ypcm 



(Bll) 



Note, that the two first terms on righthanded side of 
Eq. IjBllfl are part of V res in Eq. (|T5|l . whereas the 
third term enters into V pr in Eq l|13|l . The electrostatic 
potential V el of point multipoles can be expanded into 
series of spherical harmonics, similarly as the global CF 
potential in Eq. The corresponding contribution to 
the parameters have the following form>i& 



B 



kq 



{2k + 2n + l)l 
2™(2fc)! 



1/2 



q \x -q-fj, J Rk+n+l 



1 



cS^HRt/Rt), (B12) 



r 



where (r fe ) is a mean value of r k for the 5/ orbital. M^"" 1 
is the 2n-pole electric momentum induced on ion t, where 
t runs over the all ions in the infinite net and n runs 

(n) 

over the components of MJ . The state of electrostatic 



(n) 

equilibrium between the Mj moment and the remaining 
multipole moments of the crystal lattice is represented by 
the following equations determining the electric multipole 



moments M 



O).10.14 



M i n) = E E (~ 1 ) n+1 a * n) I(2ri) v 

t> p=0,l,2 



M (f). V W- 



(B13) 



where l( 2n ) is the diagonal unit tensor of rank 2n. oq with 
is the 2n-polc polarizability. 



APPENDIX C: 

Temperature dependence of the paramagnetic suscep- 
tibility is given by the Van Vleck formula^ 

X*(T) = ^l^ (/3a7 ct + 2& 7 , Q ) exp(-/?£ 7 ) (CI) 



a 7," = E \(i\ La +aS a \i) 



\(<y\L a + g S a \<y')\ 

E-y' E-y 



Z = ^cxp(-/?S 7) 



(C2) 



(C3) 



(C4) 
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f3 = 1/kT, a = x,y, z, index 7 numbers the eigenstates, gyromagnctic ratio of the electron spin. 
EL's denote their energies, L a and S a are the a compo- 
nent of the total orbital and spin operators and g is the 
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